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AN APPROXIMATE THEORY 
By Bernard Maggin and RoLert E. Shanks 


SUMMARY 


An experimental investigation was made to deteimine the effects of 
various design parameters on the lateral-stability characteristics of a 
glider towed by a single towline. The investigation showed that it is 
possible to obtain inherent lateral stability with a single towline system. 

An approximate theoretical analysis was also made and the results cf 
calculations made by use of this analysis were compared with the model 
flight-test results. Although the theoretical damping results are too 
conservative to be of much practical value, the existence of divergences 
and the periods of the lateral oscillations are predicted with fair 
accuracy . 


INTRODUCTION 


The Army and Navy have long been interested in towed gliders as a 
means of transporting men, material, aerial targets or guided missiles. 
One of the major problems connected with the use of towed gliders has 
been that of obtaining lateral stability of the glider on tow. Once 
trim conditions had been established, an inherently stable glider towline 
system would require no pilot attention and thus glider-pilot fatigue 
would be reduced on long flights or under conditions of poor visibility. 
In addition, it would make feasible some glider applications that are 
now impractical. 

In order to obtain lateral stability, various automatic and semi- 
automatic devices have been used but these devices are limited in their 
application by their complexity and maintenance problems. Glider- 
position Indicators, as a visual aid to the pilot during blind-flying 
conditions, have proved unsatisfactory . Various systems other than a 
single line have been proposed. One such system consisting of twin 
parallel towlines was studied in a theoretical and experimental investi- 
gation in the Langley free- flight tunnel (references 1 and 2) . In 
general, however, it appears that the most satisfacto 2 ?y solution to the 
towed-glider problem would be an inherently stable single -towline system. 
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Although a considerahle amoimt of experimental and theoretical work has 
"been conducted in this country and hy the British, to date no satisfactory 
theory predicting the lateral-^stahility characteristics of a glider on a 
single towline has Been reported. 

The results of an experimental investigation to determine the effect 
of varying the effective dihedral, directional stahllity, relative density, 
towline attachment point, and towline length are given in the present 
report. In addition, a simplified theoretical method developed to deter- 
mine the lateral-stahility characteristics of a glider towed hy a single 
towline is presented in the appendix and the resiilts of calculations made 
hy this theory are compared with the experimental data. 


The forces and moments are referred to the stahility axes (see fig. l) 
unless otherwise stated. These axes are defined as an orthogonal system 
having its origin at the center of gravity and in which the Z-axis is in 
the plane of symmetry and perpendicular to the relative wind, the X-axis 
is in the plane of symmetry and perpendicular to the Z-axis, and the 
Y-axis is perpendicular to the plane of symmetry. The coefficients pre- 
sented here refer to the glider except where otherwise noted. 

W weight, pounds 

m mass , slugs 

h wing span at zero dihedral angle, feet 

D drag, pounds j differential operator (d/ds) 

L lift, pounds^ rolling moment about the X stahility axis. 


SYMBOLS 


foot-pounds 


M 


pitching moment about the Y stahility axis, foot-pounds 


C. 


L 


yawing moment about the Z stahility axis, foot-pounds 
lift coefficient (Lift/qS) 


C 


D 


drag coefficient (Drag/qS) 



C. 



'Y 


lateral- force coefficient (Y/qS) 
weight coefficient (w/qS) 
rolling-moment coefficient (L/qSh) 
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pitching-moment coefficient (M/qSc) 

Cn yawing-moment coefficient (N/qSh) 

^ angle of yaw, radians 

\|r* angle between X-axis and projection of the towline on 

the X-Y plane 

3 angle of sideslip, radians 

(ji angle of hank, radians 

9^' angle between Z-axis and projection of the towline on 

the Y-Z plane 

p rolling angular velocity, radians per second 

r yawing angular velocity, radians per second 

( angle between towline and the relative wind, degrees 

a angle of attack, measured from the top of the boom to the 

relative wind, degrees 

angle of attack, measured from the angle of zero lift, degrees 
r wing dihedral angle (positive wing tips up) , degrees 

Lp rate of change of rolling moment with angle of sideslip (SL/^P) 

Lp rate of change of rolling moment with rolling velocity (Sl/^p) 

rate of change of rolling moment with yawing velocity (^L/Sr) 

Np rate of change of yawing moment with angle of sideslip (Sn/^P) 

c mean aerodynamic chord of wing, feet 

S wing area, square feet 

vertical- tail area, square feet 

t time, seconds 

p mass density of air, slugs per cubic foot 

q dynamic pressure, pounds per square foot ( lp 
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gravitational acceleration , ft/ sec^ 


relative-density factor 


(; 


pSt 


m 




radius of gyration about X-axis, feet 
radius of gyration about Z-axls, feet 

radius of gyration about X-axis, spans 
radius of gyration about Z-axls, spans 

radius of gyration about principal longitudinal axis, spans 
radius of gyration about principal normal axis, spans 


angle of attack of principal longitudinal axis of airplane^ 
positive when forward end of major principal axis is above 
X-axis, see figure 1 

distance along X-axis from center of gravity of glider to 
the towline attachment point, spans 

distance along Z-axls from center of gravity of glider to the 
towline attachment point, spans 

sidewise movement of center of gravity along Y-axis, feet, 
see figure 8 

sidewise movement of center of gravity along Y-axis, spans 

towline length, spans 

true airspeed, feef per second 

velocity along Y-axis, feet per second 

towline tension, pounds 


product-of -inertia factor, spans' 


,2 



lateral force, pounds 
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rate of change of lateral- force coefficient -with angle of 

/aCyX 

sideslip, per radian ( 

rate of change of rolling- moment coefficient with angle of 
sideslip, per radian (SCj/Sp) 

rate of change of rolling-moment coefficient with rolling 
angular velocity factor 


velocity factor 



rate of change of rolling-moment coefficient with yawing 

2.1 > 


rate of change of yawing-moment coefficient with angle of 

/Sc. ' 


n 


sideslip, per radian I 


rate of change of yawing-moment coefficient with rolling 
angular velocity factor 





C rate of change of yawing-moment coefficient with yawing 

angular velocity factor ( 

\ 2V> 

A, B, C, E, F, G, H, I, J, K coefficients of the stability equation 
P 

h 

2 

id 


period of oscillation, seconds 

time required for a motion to damp to one-half amplitude 


imaginary portion of complex root 
c real portion of complex root or a real root 

Towline terms : 

Y rate of change of the lateral force with sidewise displacement 

y (dY/Sy) 


Yj 


rate of change of lateral force with angle of yaw (Sy/S\|t) 
rate of change of lateral force with angle of bank (dY 
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rate of change of rolling moment with sidewise displacement 

(^L/^y) 

rate of change of rolling moment with angle of yaw 

OL/Si) 

rate of change of rolling moment with angle of hank 

rate of change of yawing moment with sidewise displacement 

(^N/Sy)- 

rate of change of yawing moment with angle of yaw (Sn/^tIt) 
rate of change of yawing moment with angle of bank (dN/B0) 


rate of change of lateral foixe caused hy towline tension 
with sidewise displacement (^Ty/^y*) 

rate of change of lateral force caused hy towline tension 
with angle of yaw 

rate of change of lateral force caused hy towline tension 
with angle of hank 

rate of change of yawing moment caused hy towline tension 
with sidewise displacement 

rate of change of yawing moment caused hy towline tension 
with angle of yaw ) 

rate of change of yawing moment caused hy towline tension 
with angle of hank 

rate of change of rolling moment caused hy towline tension 
with sidewise displacement 

rate of change of rolling moment caused hy towline tension 
with angle of yaw 

rate of change of rolling moment caused hy towline tension 
with angle of roll (bT^/d0) 
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TESTS 

APPAEATUS 


All the tow tests of the model were conducted In the Langley free- 
flight tunnel, a complete description of which is given in reference 3* 
Photographs of the model mounted on a stand and on tow in the tunnel are 
presented in figure 2 . 


A sketch of the model used in the tests is presented in figure 3* 
The model consisted of a wooden boom upon which the wing and stabilizing 
surfaces were mounted. A Ehode St. Genese 35 airfoil section was used 
in accordance with free-flight- tunnel practice of using airfoil sections 
that give maximum lift coefficients in low- scale tests approximately 
equal to those obtained in fxill-scale tests of conventional airfoil 


sections. The stabilizing surfaces were constructed of inch sheet 


balsa . 


The wings were mounted so that a range of geometric dihedral angles 
between and 15 ° could be obtained, and the model was arreinged to 
allow for mounting of any one of three vertical tails, 5 percent, 

7.5 percent, and 10 percent of the wing area. (See fig. 3«) 

The glider model was equipped with conventional control siirfaces 
actuated by a pilot through standard free-flight-tunnel control 
mechanisms. A complete description of the flight models and flight 
technique used in free-flight-tunnel testing is given in reference 3. 

The model photograph and sketch in figures 2 and 3, respectively, 
show the special tow bar used on the model. The tow bar consisted of a 
horizontal bar mounted on two vertical bars, one located in the nose of 
the model, the other at about the 60 -percent station of the root chord 
of the wing. The vertical position of the horizontal bar could be 
varied and the towline attachment could be made anywhere along the 
horizontal bar. The tow bar was used for convenience, but in practice ‘ 
a trifurcated bridle system, which gives, in effect, a fixed attachment 
point at the apex of the glider bridle lines could be used, provided 
that the attachment lines remain in tension. Figure 4 shows some of 
the more commonly used single -towline attachment systems. From a study 
of these systems it can be seen that the tow-bar arrangement used in 
the tests can represent any of these common systems insofar as the towline 
attachment point is concerned if the longitudinal trim is assumed not 
to change (as is approximately the case) during a disturbance. 

The relative density p of the model was adjusted without changing 
the moments of inertia by adding weights at the model center of gravity. 
Increasing the value of p in this manner decreased the radii of 
gyration as shown in table I. 
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METHOD 


In the tests the airspeed was held constant and the trim angle of 
the elevator was adjusted to obtain approximately the desired angle of 
attack and towline angle. The model was controlled by the pilot who, 
in the case of stable tow configurations, supplied disturbances to start 
an oscillation. In the case of unstable tow configurations, the pilot 
steadied the model so that the oscillation could develop from the steady 
state and then could be stopped when sufficient records of the motion 
had been obtained. During the oscillation the controls were fixed. 

The motions of the model for the various test conditions were 
recorded by two motion-picture cameras, one mounted on top of the tunnel 
directly above the glider model, the other at the rear of the tunnel 
directly behind the model. In some instances the period of the lateral 
motion was measured with a stop watch. For each tow condition the tow- 
line angle and angle of attack were measured visually with a protractor 
mounted at the side of the tunnel. 

From the mot ion- picture records, plots were made of the sidewise 
displacement, angle of bank, and angle of yaw against time for the con- 
ditions tested. From these records the period and time to damp to one- 
half amplitude of the lateral oscillations were obtained. Some repre- 
sentative flight records are presented in figures 5(a) and 5(b) . 


SCOPE 


A list of the conditions covered in the tests is presented in 
table II. It will be noted that the tow attachment points considered 
are above and forward of the glider center of gravity and the glider is 
below and behind the tug. These glider tug configurations were used 
because past towing experience indicated that they would be best from 
the standpoint of obtaining lateral stability on tow. This table shows 
the effects on lateral stability of varying the effective dihedral, 
directional stability, relative density, towline attachment point, and 
towline length. In order to determine the effects of these parameters. 
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they were, varied one at a time from a basic condition for which the 
values of these parameters were: 

W = 0.1)-65 lb u = 2.i+ 

I = 4 

b 

^ = 0.558 

b 

^ = 0.225 

b 

r=io° c, =-0.15 

S+ 

-^ = T ‘5 percent Cj, = O.06 

S 

A few tests were also made to detenaine the effect of varying the effective 
dihedral for zero-length towline. For these tests the model had three 
degrees of freedom about the tow attachment point. The values of the other 
parameters were the same as for the basic condition given above. 

Test results are available only for towline lengths up to 4 spans 
because of the size of the tunnel test section. Analysis of the towline 
terms and calculations extending the range of towline length indicate 
that changing the length of the towline up to about 10 spans affects the 
lateral- stability characteristics appreciably, but further increase has 
only slight effect. Therefore, for towline lengths of 10 spans and 
greater, the trends which are found for the remaining parameters at the 
shorter towline length can be expected to prevail. 


THEOEETICAL CALCULATIONS 


In addition to the experimental investigation an attempt was made 
to develop an approximate theoretical method for predicting the lateral- 
stability characteristics of the glider on' a single towline. By use of 
this theoretical method, which is given in the appendix, calculations 
were made for various combinations of towline length, effective dihedral 
parameter 


C^p, directional stability parameter 


■^np, relative -density 


factor u, and 'vertical and horizontal towline attachment position. A 
complete listing of the conditons for which the calculations were made 
is given in tables I and II . Inasmuch as the product of inertia was 
believed to be relatively unimportant for the low-relative-density 
condition of most of the model tests, the product- of- Inertia factor 

was assumed to be zero for most of the calculations . 
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The parameters used In the calculations were obtained either from 
tests or from C6ilculations for the specific model condition. The source 
of each of these parameters is Indicated in table I. Except in the 
relative-density derivatives, values of the derivatives for an angle of 
attack a of 2° and a towline angle € of 25 ° were used in the 
calculations . 


RESULTS AND DISCUSSION 


The results of the investigation are presented in table II and in 
figures 6 and 7 in which lateral oscillations are presented in terms of 
the period and the reciprocal of the time to damp to one-half amplitude 
for each particular condition. The reciprocal of the time to damp to 
one-half amplitude is used to evaluate the damping because this value 
is a direct measure of the degree of stability. In general, the calcu- 
lations and tests were made at the same conditions except where the 
dynamic characteristics of the model or physical limitations of the test 
setup precluded testing. 

The lateral- stability chsiracteristics determined by the theory pre- 
sented in the appendix indicate that the lateral motions of the glider 
model generally consist of two highly damped aperiodic modes and two 
periodic modes, one usually having a period of approximately three times 
the other. From the results of the calculations and tests it appears 
that the long-period mode is usually the predominant or more lightly 
damped motion and consists of a combined rolling and yawing motion 
slmilEir to that of the conventional Dutch roll. 

Wherever possible the lateral- stability characteristics of both 
periodic modes were obtained from flight records similar to those of 
figure 5* In general, the characteristics of the predominant long- 
period mode were obtained from the plot of sidewise displacement against 
time (fig. 5(a))* The short-period mode appeared primarily as a 
yawing oscillation and its characte id. sties were obtained from plots of 
yaw angle against time (fig. 5(b)). In most cases it was difficult 
to ascertain the period and damping of the short-period oscillation 
because it was masked by the less heavily damped long-period oscillation. 

Althou^ no systematic investigation was conducted to deteimlne the 
effect of towline angle on the lateral- stability characteristics of the 
glider model, analysis of the towline forces Indicated its Importance. 

Some exploratory tests were made to verify this analysis and it was found 
that, in general, for any attachment point in the range tested, increasing 
the towline aiigle was stabilizing. During the tests an attempt was made 
to hold the towline angle and angle of atte.ck to one set of values 
( € = 25°^ a = 2°) . Although these values of € and a were not always 
exactly obtained, it is believed that the slight variations of these 
parameters do not invalidate the correlation of theoretical and test 
results . 
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COEEELATION OF CALCULATED AND EXPERIMENTAL RISULTS 


The results of the experiments and analysis which are presented in 
table II and figures 6 and 7 indicate that, althou^ the theory predicts 
the periods of the lateral oscillations fairly well, it does not predict 
the damping with sufficient accuracy to be of practical use • 

In order to determine whether the omission of the product- of- inertia 
terms, recently found to be important in some cases (reference 4), was 
responsible for the poor quantitative agreement of the calciilated and 
experimental damping results, additional calculations were made for the 
basic condition and for the relative-denslty-of-10 condition with these 
terms included in the equations • These calculations showed that inclusion 
of the product- of- Inertia terms had virtually no effect on the results 
for the basic condition (u = 2.4) but did change appreciably the damping 
results of the higher relative-density condition (U = 10 .0). Therefore, 
with the exception of the relative-density variation, product- of- inertia 
effects were ignored for these tests and calculations . 

In view of the fact that only first-order effects were considered 
in developing the theory, the discrepancies between the theoretical and 
experimental results may in part be associated with the relative importance 
of some terms which were considered negligible to the first approximaticn . 


Effect of C^ 
^3 


of 


Four- span towline le 
given in figure 



The experimental data for the range 
show the periods of both the long-period 


and the short-period modes to be fairly constant. Although the short- 
period calculations are in good agreement with the tests, the long- 
period results Indicate increasing length of period with decreasing 
effective dihedral. 


% 


The experimental damping results Indicate that the long- period motion 
Is the predominant mode since the short-period motion was always very 
heavily damped. The damping of the long- period oscillation was found to 
decrease with decrease in ^ value of -O .0458 

(table II) any disturbance results in a rapid divergence in roll which 
obscures the characteristics of the oscillatory modes . The calculated 
results do not agree with the teste except in the prediction of the 
divergence . 


Zero towline length .- The test results of figure 7 show that for 
the zero towline length the short-period oscillation is the predominant 
motion for all positive values of "C^p. The length of the period, 


however, increases as the value of 



decreases indicating a 
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transition of the predominant motion from the short-period to a longer- 
period mode at low and negative effective dihedrals . The calculations 
show the short-period motion to he the predominant mode only at the higher 
values of that the long period is the more lightly damped at 

lower effective dihedrals . These results also clearly Indicate that the 
long- period mode is replaced hy a divergence as the value of -C. is 

P 

increased from 0 to -0.0458. 


Effect of C 


n, 


P 


The test results of figure 6(h) show little change in the period of 
the long- period oscillation for the range of C tested, whereas the 

period of the short-period oscillation increases gradually with reduction 


in C 


n, 


P 


At zero directional stability the short-period motion which was 


hardly noticeable at higher values of became evident as a steady 

P 

large -amplitude short-period yawing oscillation. It was assumed that the 
motion was unstable at small amplitudes and built up to a steady oscilla- 
tion of large amplitude . The rate of increase of amplitude could not be 
measured, however, because the oscillation was well developed by the time 
the trim conditions had been attained and the short-period yawing motion 
could not be stopped by the 'pilot to permit study of the motion at small 
amplitudes . The calculations show a marked decrease in the damping of 
the short-period mode with decrease in and at C =0 indicate 


P 




the unstable oscillation which was found experimentally. 


With increase in the directional stability the damping of the long- 
period oscillation was found to increase slightly. The calculated results 
are not in agreement with this trend. 


Effect of \i 

An indication of the effect of the relative -density factor [i on 
the lateral- stability characteristics is given in figure 6(c) . Because 
of test limitations imposed by the model, the angle of attack as well as 
airspeed had to be increased as the value of }i was increased. Table I 
gives the corresponding variation in the other parameters used for the 
calculations. Compared to the change in magnitude of the values of u, 
the changes in the other parameters were relatively small, however, and 
the principal effect given in figure 6(c) is believed to be that of 
variation of \i . 

Although the period of the long- period oscillation increases slightly 
with increase in relative density, the long-period calculated results show 
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a very slight decrease in length of the period. The inclusion of the 
product- of- inertia tem In the calculations for a value of p. = 10 had 
practically no effect on the period of either the long-period or the 
short-period oscillations- The experimental results show a general 
Increase in damping of the long-period mode as the value of Is 
increased from 2.4 to 10, hut the calculated results presented in 
figure 6(c) indicated that increasing [i from an average present-day 
value of 2.4 to 10.0 results in a change from stability to gradually 
increasing instahlllty for the long-period oscillation- Including the 
product- of- inertia teims in the calculations (table II) Improves the 
agreement between damping results of the tests and calculations at a 
value of [1 of 10 but produces virtually no change in the results at 
a value of.\i of 2 -4 . Since the tests gave no indication of short- 
period instability at any value of relative density, it is apparent that 
including' the product- of- inertia terms in the calculations made the 
agreement between the tests and calculations poorer for the short- 
period mode - 


Effect of Tow Attachment Point 

The effect of varying the horizontal location of the tow attach- 
ment point was not determined experimentally because the model could 
not be trimmed longitudinally for values of x below that of the basic 
condition- However, calculations were made for a range of values of x 
from 0 to 0.558 spans, and these results are presented in figure 6(d) 
with the one experimental point available . The results show a decrease 
in the instability^ of the long- period motion and a decrease in damping 
of the short-period motion with increase in x- With increase in x, 
the period of the long-period oscillation decreases at a decreasing 
rate and the length of the short-period oscillation is practically 
constant- The effect of the vertical location of the the tow attach- 
ment point is given in figure 6(e). Increasing the z distance was 
found to have little effect on the length of the period of either mode 
but reduced the length of the longer- period mode slightly. The test 
results show that the greatest damping was found at a value of z of O.II 7 
and that increasing or decreasing the horizontal distance reduced the 
damping of the long-period mode. The calculations indicate, however, 
that increasing z improves the stability of the long- period mode but 
reduces the short-period damping. 


Effect of Towline Length 

Both the theoretical and experimental results presented in figure 6(f) 
indicate that increasing the towline length from 1 to 4 span lengths 
increases the damping of the long-period oscillation. The test results, 
however, show somewhat more damping than is predicted by the theory. 
Calculations were also made for towline lengths of 10 and 100 spans 
although tests were not possible for these lengths. The long-period 
damping continues to increase (table II) but at a much more gradual rate 
than was found for the range tested. 
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The calculated results show that increasing the towline length from 
1 to 100 span lengths has no appreciable effect on the damping of the 
short-period oscillation. 

An increase in the towline length causes a gradual increase in the 
period of the long-period oscillation hut causes no change in the period 
of the short- period oscillation. These results are shown hy both the 
calculated and experimental results which are in good agreement. 


CONCLUDING EEMAEKS 


The results of the experimental investigation and theoretical 
analysis to determine the lateral- stability characteristics of a glider 
towed by a single towline may be summarized as follows: 

1. The investigation showed that it is possible to obtain inherent 
lateral stability with a single-towline-glider arrangement. 

2. The stability theory presented in the present paper does not predict 
the stability of the glider model with sufficient accuracy to be of much 
practical value . The calculated and experimental results are in fairlyr 
good agreement on the period of the lateral oscillations and the -prediction 
of divergences, but the calculated damping of the lateral oscillation is 
generally considerably less than the measured damping. The theory is 
presented, however, as a basis for further study and the experimental 
results should be useful for correlation with any further theoretical work 
which may be performed. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 
THEOEETICAL METHOD 


The method of calculating the lateral atahlllty of a glider towed 
hy a single towline consisted of setting up the stability equations with 
respect to the stability soces of~ the glider (see fig. l), obtaining the 
aerodynamic stability derivatives for the glider and the derivatives 
caused by the towliiie forces, and solving for the period and time to 
damp of the lateral oscillatory modes of the glider. 

Assumptions •- To simplify and facilitate the handling of the theo- 
retical analysis, the following assumptions were made: 

1. The glider and tug were assumed to be in level flight and the 
tug to be in steady flight. 

2. The basic aerodynamic peirameters were assumed to vary independently 
of each other. 

3- Lateral- stability characteristics were assumed to be independent 
of the longitudinal-stability characteristics . 

14^. The angles of deviation of the glider from the steady-flight 
condition were assumed to be small. Towline angles with respect to the 
relative wind and towline tension were assiimed not to change during a 
disturbance . 

5. The towline was assumed to be weightless and straight and to ha-ve 
no effect on the stability of the glider other than that caused by direct 
forces and moments resulting from towline tension. 

Full-scale tests and tests in the Langley free-flight tunnel have 
indicated that the assumption of a straight towline is valid for moderate 
towline lengths . British investigations show that glider towlines can 
be considered strai^t for lengths up to about h or 5 spans for any 
towline angle and for greater lengths for towline angles between about 20° 
to 30°. 

Equations of motion. - Inasmuch as a single towline allows the glider 
three degrees of lateral freedom, the basic lateral- stability equations 
were used, modified only to include the derivatives defining towline 
restraints. The lateral equations of motion following a disturbance may 
be written as follows: 

m^ + mVr = mg0 + 3Yp + towline derivatives (l) 

Tok'^ ^ ~ ^ ^^p towline derivatives 


( 2 ) 
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m]r 2 ^ ~ + + PLo + towline derivatives ( 3 ) 

dt 2V P dt 2V ^ P 

If the effect of glide-path angles other than zero is considered, 
the normal glide-path terms will appear in the basic lateral-stability 
portion of the preceding equations • 

Towline derivatives .- In straight steady flight the towline con- 
tributes only a pitching moment and does not affect the lateral stability. 
Sidewise displacement of the glider along the Y-axis, displacements in 
yaw about the Z-axis, or displacements in roll about the X-axis result 
in forces and moments about these axes that are functions of towline 
tension.. The towline derivatives expressing the resultant relationship 
may be obtained to the first order as follows from figure 8 by consid- 
eration of the forces and moments with respect to the stability axes 
and the glider center of gravity: 

(a) Lateral-force derivatives 

1 = -1 
T I 


or 


and 


Y = 

y 


T 

I 


Y _ _ X sin i 
T I 


(cos G sin \lf') 


but for relatively long towlines, however, ijr* = then 


where sin ''I'’ ' = and 


— = sin (i - sin e sin (i 
T I 


Also for long towlines, = (jS or 


V = -T(f + 0 




b) 


( 6 ) 


where sin (j( = (^ . 

'(b) Yawing- moment derivatives: 



(7) 


NACA EM No. L8H23 


17 


N 




= -Tx^^ ^ cos ^ 


N^ = “Tx^ + sin € 
(c) Rolling-moment derivatives: 

^ = 


I 


L. = -Tz 




= -Tz 


(! 

(! 


+ cos € 


+ sin €' 




') 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


where T = 


Draf^ 
cos € 


and 


€ = tan' 




For zero towline length, the tow attachment point is fixed relative 
to the tug and the towline tension always acts in the vertical plane of 
the relative wind- Because of these limitations the glider is permitted 
only three degrees of freedom and the towline derivatives Yy, Ny, and 

normally resulting from sidewise displacement of the tow attachment 
point do not exist. The remaining derivatives resulting from yaw and 
roll about the tow attachment point become: 


(a) Lateral-force derivatives : 


= -T cos e 
Yw = -T sin 6 


(b) Yawing- moment derivatives: 


N, = -Tx cos € 
Nw = -Tx sin e 


(c) Rolling-moment derivatives: 


L, = -Tz cos G 


(5a) 

(6a) 

(8a) 

(9a) 

(11a) 


= -Tz sin 6 


(12a) 
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Equations of motion, nondlJiienslonal form.- Substituting the towllne 
derivatives In eq,uatlons (l), (2), and (3) yields the following eq.uatlonS : 


ng + mVr - mg0 - pYp - yYy - tY,|r - 0Y0 = 0 (13) 

mkz^ g - mixz || - ||Np - ~ 3N3 - yNy - - 0N0 = 0 (ll^) 

ft ■ ~ df ■ - PLp - yLy - - 0L0 = 0 (15) 

For convenience, the equations can be converted to nondlmenslonal 
form, by using the following relationships and the operator D : 


r' = I 


= 


HI 

pSb 


— TT 


K. . vS 


y = V J (P + ilf) dt = 

^ = V 

dt 

^ = p 

dt P 

dt 


dt 


= r 


P + t 


- = P 

V 


D = f- 
ds 


where 
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and 


ds = J dt 

2 2Y 

54f _ rb 
2 2T 

Then substituting these relationships and dividing equation (l3) 
by ^SV^b ■ yields the nondimensional equations : 

(2W)3 - CYpI^ - Ty^.o) p + (i«.d 2 - 2Ty^D - aTy^j) § 

* (-20„I) - 2Ty^D) §.0 ( 16 ) 


-Cngl^ - Tny.D) P + ^IwKz^dS - - 2Tn^D - 2T„y,j g 


+ i'-MjzDS - C^d 2 - 2Ty>) = 0 


(17) 


- TZy.D) P + - Cj^d 2 - 


2Y 


+ - CZpD? - 2 Tj0d) ^ = 0 


(18) 


vhere the nondimensional forms of the towline derivatives are 


Ot 

= ~T^ ^eing 
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. = -xT 






T = -xT 
n 




T = -vT 


T, = -zT. 


^y' 


= 


■zT„ 


The determinant of the left side of these equations yields the 
glider-stahility equation of the foim: 

AD^ + BD® + Cd'^ + ED^ + Fd5 + gD^ + HD^ + ID^ + JB + K = 0 

where : 


(19) 


A = 32^3k^2k^2 _ 32h3Ky^2 


■xz 


B = 


-8u ^ 2 . ^XZ*^np 

- 8n^Ey 




C = 


2o ^ 1. 2 


9( ^ ^ ^ •'y' ■ ^ Zp 

2 „ 2 , 

+ '+|ajv-5f uy v-n + lt>n 
3 "p ^ ^3 r + - -p 




*■ ^XZ ‘ ^Xz'^i^ 
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E 


' - "HiCn T. ‘ 'wCj t_ + T, 


"r V 
2 


■ -p ^ ^3 


“«^Z * “H^S: Vj,' ■ '^n T 


P h 


^ \\\ ■ '^w * 


•» “p 


F = 


SyiK^CiC^ - 

- 2 H C T, + 2C C T, - 2n C, T + 2 Cy„C, T„ 

Yp ^ ^ H ^P ^p ^ ^r (jf 

- Cn C^ T + c Cj T - 8 nC T- + 8 nC T + 2 Cj^ C, T 

^ p ^y' ^P r ^y' ^p V ^P 3 \ 

8^Ky. Cj^pTy^, - 8uK^Ty^,C2p - 8MKx2;CvTn^, + 8[iKj2pY , 

- "C2pC^C^ - 2C C2^C^ 


^ ■ ^^ZpVn^‘ + ^%VZyi ■ ^Cx^C^Tn^, - 


H = 


+ 20 C T, + 2 C C, T 

"'p"pV “s p^y' 


r y 


P P y 


2 C, C_ T 
Z n y 

P P ^y 


I = j = K = 0 
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The period and time to damp to one-half amplitude of the oscillatory 
and aperiodic modes may he determined from the roots of equation 19 hy 
the following relationships 

p _ 2rt h 

d V 

T = h 

1 c V 

2 

or 

m _ 0 .0693 h 

1 ■ c V 

2 

where d is the coefficient of the imaginary portion of the complex 
roots and c is either the real portion of the complex roots for the 
osclllatoiy mode or the real root defining the aperiodic modes. 
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TABLE I.- PHYSICAL AND AERODYNAMIC CHARACTERISTICS 
OF THE SINGLE TaCLINE GLIDER MODEL 


Wing area, square feet 1,02 

Wing span, feet 2.50 

Wing aspect ratio , 6.10 

Wing mean aerodynamic chord, feet 0 . 4 l 6 

Model cente 3 X)f— gravity location, percent mean chord 0,390 


Mass characteristics 
(All values obtained by measurement) 

Weight 

W 

Relat ive-dens i ty 
factor 


% 


""p 

^xz 

^211 

1^50 

720 

900 

2,4 

5.0 

8.0 
10.0 

0.419 

.364 

.289 

.259 

0.462 

0.606 

.543 

.431 

.388 

0.596 

-0 .0240 




.259 

.388 

- .0143 


Aerodynamic characteristics 


% 

Cl 

Qd 

cir 

‘^np 

% 


(deg) 

Cb) 

0>) 

(c) 

(cf 

(c? 

^ 2.4 

10 

0.57 

0.110 

0.161 

-0.0272 

-0.49 

5.0 

12 

.75 

.150 

.221 

-.0374 

-.49 

8.0 

Ik 

.86 

.168 

.282 

-.0476 

-.49 

10.0 

15 

.92 

.175 

.322 

-.0544 

-.49 


Vertical 

tail 

Area 

(percent S) 

Cnp 

(t) 

Cyp 

(■b) 

Cnr 

(b) 


0 

0 

-0.257^ 

- 0.020 

1 

5.0 

.0343 

-.4061 

-.040 

a 2 

7.5 

.0572 

-.4462 

-.060 

3 

10.0 

.1087 

-.5205 

-.082 


Dihedral angle 
(deg) 

(b) 

-5 

0.0458 

0 

0 

5 

—.0688 

°10 

-.1375 

15 

-.2000 


®'Basic condition 
^Values obtained 
^Values obtained 


from which the parameters were varied, 
from force tests in the Langley free— flight tunnel, 
from reference 5 . 


TABLE n.- CORRELATION OF THE THEORETICAL AND EXPERIMENTAL RESULTS OF THE LATERAL 
STABILITY INVESTIGATION OF A GLIDER TOWED BY A SINGLE TOWLINE 



®Neutrally stable from visual observation. 


ro 

VJl 
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Y 




direction 




Figure 1.- The stability system of axes. Arrows indicate positive directions 

of moments, forces, and velocities. 
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(b) In flight. 

Figure 2.- The glider model used in the tow tests. 
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Figure 3. - Sketch of the glider model. 
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Figure 4.- Glider tow systems employing a single towline. 


Sjdei^ise displgcemenf^ fl 
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A' 




rs 

1 

A 


1 \ 



V 




\\ 



1- fspan 
1 1 ^ 


V, 






(a) Sidewise displacement. 





Recorded mohon 
Amplitude 
em/ebpe 


(b) Yaw angle. 

Figure 5.- Typical records of lateral motions as obtained from motion- 
picture records of tow tests of the glider model. 
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PERIOD 
OF MODE 

Long 
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Measured 
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TESTS 
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observation 


Basic condition from 
which the parameters 
were varied 
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Product of inertia: 
neglected included 
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A 


o 

0) 

la 


xi 

o 
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01 
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H 
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o 

o» 

to 


o 

•H 

Q> 

PU 


V. 

SI 


o 

stab] 


O 

Unstable 

1 



0 -.08 -.16 .24 
Effective dihedral 
parameter 9 C^o, per 
radian ^ 




ment point location, x, 
spans 

(d) 




O 



N 






^'st 

»bla° 

’7777^ 

/■ 

/ 

1_ 

Unstable 

L_] 


0 .04 .08 .12 


Directional stability 
parameter, Cuq, per 
radian ^ 


(b) 


0 

o 

) 0 
A 




B- 

— 

-A 





tA 

Ostaf 


// / / 

" l/nsta^e 


0 4 8 12 

Relative density 
factor 


(c) 
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o 





— 

□ 

-H 





o s 

o ^ 

stable 


//y/At 


7/// 


Unstable 


0 .10 .20 .30 
Vertical tow attach- 
nont point location, 
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(©) 





ern' 

: 

- — 


n 

r 


o< 


— — 
Hble 


y//// 

Qnj 

'7777 

liable 

1 
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L e 

^ u 


Towline length 
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Figure 6.- Correlation of the period and the reciprocal of the time to damp 
to one -half amplitude of the lateral oscillations of the glider model as 
obtained from tests and theory. (See table II for complete conditions of 
each part. ) 
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Long period 
Short period 


o 

Q> 


PU, 




o 





Dihedral effectiveness 
parameter I > per radian 


Stable 


Figure 7.- Correlation of the period and the reciprocal of the time to damp 
to one -half amplitude of the lateral oscillations of the glider model at 
zero towline length as obtained from theory and tests for various values 
of ~Ci (For other conditions see table II,) 
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(a) Steady flight condition; « = 0,y = 0. 





Tunnel screen 


(c) Yaw. 



Figure 8.- Sketch showing the relationships used in obtaining the towline 

derivatives. 



